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Abstract 

We calculate final-state interaction corrections to the energy-angular distribution of in 
semi-leptonic top quark decay, where the parent top quark is produced via e+e" — > it near 
threshold. These are the corrections due to gluon exchange between t and b {i and b) and 
between b and b. Combining with previously known other corrections, we explicitly write 
down the energy-angular distribution including the full 0{as) = 0{P) corrections near 
tt threshold. Numerical analyses of the final-state interaction corrections are given. We 
find that they deform the distribution typically at the 10% level. We also find that 
all qualitative features of the numerical results can be understood from intuitive pictures. 
The mechanisms of various effects of the final-state interactions are elucidated. Finally we 
define an observable which is proper to the decay process of the top quark (dependent only 
on dTf-^jji+j^/dEidQi of a free polarized top quark) near tt threshold. Such a quantity will 
be useful in extracting the decay properties of the top quark using the highly polarized 
top quark samples. 
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1 Introduction 



A future e^e~' linear collider operating at energies around the ti threshold will be one of the 
ideal testing grounds for unraveling the properties of the top quark. So far there have been a 
number of studies of the cross section for top-quark pair production near the tt threshold, both 
theoretical and experimental |]ll-|ll9|, in which it has been recognized that this kinematical 



region is rich in physics and is also apt for extracting various physical parameters efficiently, 
e.g. mt, as, Tt, rriH, QtH, etc. 

While most of the previous analyses were solely concerned with the production process 
of the top quark, one may also analyze the decay process in detail and extract some important 
physics information. Especially the fact that t and i are produced highly polarized in the 
threshold region is potentially quite advantageous for studying the electroweak properties of 
the top quark through its decay. Detailed investigations of the decay of free polarized top 
quarks have already been available including the full 0{as) corrections [jl6|, Close to 

threshold, however, these precise analyses do not apply directly because of the existence of 
corrections unique to this region, which connect the production and decay processes of the top 
quark. Specifically, these are the corrections due to gluon exchange between t and h {t and h) 
or between h and h. 

This type of corrections arises when the particles produced decay quickly into many 
particles, and are referred to as "final-state interactions", "rescattering corrections", or "non- 
factorizable corrections" in the literature. They generally vanish in inclusive cross sections 
1^, |T^, |TT|, but modify the shape of differential distributions in a non-trivial manner [^, |T3], |T^|. 



The size of the corrections is at the 10% level in the tt threshold region, hence it is inevitable to 
incorporate their effects in precision studies of top-quark production and decay near threshold. 



The same kind of effects has recently been studied in W pair-production |2^, |2^, |29 . 

The first analysis of the top-quark decay in the threshold region was given as a part of the 
results in Ref. |19[. In that paper, the mean value {nt) of the charged lepton four- momentum 
projection on an arbitrarily chosen four-vector n in semi-leptonic top decays was proposed as 
an experimentally observable quantity sensitive to top quark polarization, and this quantity 
was calculated including the final-state interactions. Clearly, and also admittedly in that work, 
the calculation of the differential distribution of l'^ including the final-state interactions has 
been demanded. 

In this paper, we calculate these final-state interaction corrections to the differential 
energy- angular distribution of the charged leptons. We find that the corrections deform the 
distribution non-trivially at the expected level. Combining with other, previously known 
results, we write down the explicit formula for the energy-angular distribution including all 
0{as) = 0{f3) corrections near ti threshold. 

Another aim of this paper is to present physical descriptions of the final-state interactions 
which enable us to qualitatively understand the features of our numerical results. Such de- 
scriptions would be useful since the systematic calculation of final-state interactions based on 
quantum field theory is rather complicated, involving box- and pentagon-type diagrams, and 
it is not easy to make physical sense out of the obtained final expressions. To our knowledge 
such qualitative explanations have never been put forth, although corresponding theoretical 



1 



calculations and numerical studies have been partly available. 

Finally we propose an observable which is proper to the decay process of the top quark 
near threshold. Since the final-state interaction connects the production and decay processes of 
the top quark, it destroys the factorization property of the corresponding production and decay 
cross sections. In order to study the decay of the top quark in a clean environment, it would be 
useful if we could find an observable which depends only on this process {dVt^i,i+y/ dEidVLi of a 
free polarized top quark). In fact, such an observable can be constructed, which at the same 
time preserves most of the differential information of the Z+ energy- angular distribution. 

In our numerical analysis we solve the Schrodinger equation numerically in order to in- 
clude the QCD binding effects near threshold. We follow both the coordinate-space approach 
developed in Refs. 0, || and the momentum-space approach developed in Refs. [§, |1 in solving 
the equation, and compare the results. There are small differences in the numerical results 
obtained from the two approaches, reflecting the difference in the construction of the potentials 
at short distance. (The difference is formally 0{a1), of the order beyond our present scope.) 
This issue is also discussed. 

In section 2 we introduce some notations to be used in later sections. Section 3 contains 
the physical descriptions of the effects of final-state interactions. The results of the systematic 
calculation of final-state interaction corrections to the energy- angular distribution, as well 
as the complete formula for the distribution including all 0{as) = 0{f3) corrections, are given 
in section 4. Section 5 shows various numerical results and a comparison with the qualitative 
picture. We define an observable proper to the top decay process in section 6. Discussion and 
conclusion are presented in sections 7 and 8, respectively. 



2 Definitions and Conventions 

We consider longitudinally polarized e+e~ beams throughout our analyses. Pe± denotes the 
longitudinal polarization of e^, and we set 

We choose a reference coordinate system in the tt cm. frame. Three orthonormal basis vectors 
are defined as 

a„, (2.2) 



Pe- 


Pe 


X Pt, 




nil = j- 


, nN = -. 




, n_L = hn 


|Pe-| 


Pe- 


X Vt\ 





where Pe- and pj represent the e and t momentum, respectively. 

Our conventions for the fermion vector and axial-vector couplings to the Z boson are 

= 2/| -4g/sin2 0w, = 2//, (2.3) 

respectively. Certain combinations of these couplings will be useful below: 

fli = 9e?t + (^e + o.l)vtd'^ + 2qeqtVeVtd, 
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and 

c±{x) - 

Finally, we define 



a2 
as 

d 



«2 + XO-l 
0-1 + XO'2 ' 

1 04 + 

2 ai + ' 



2veaeV^ d + 2qeqtaeVtd, 
Avea^vtatd^ + 2qeqtaeatd, 
2{vl + al)vtatd^ + 2qeqtVeatd, 

1 s 
16 sin^ cos2 s - M| ' 



[2A) 



ciix) ^ 

C^n(x) = - 



2\ Ct2'^3 — CilCi4 



1 03 + xa4 



2 



2 ai + X0.2 



21/ 



-Cfb(x). 



:2.5) 



2 ' 

mf 



l + 2y 



(2.6) 



The QCD enhancement of the top-quark production cross section near threshold is in- 
corporated through the S'-wave and P-wave Green's functions, G{p,E) and F{p,E), of the 
non-relativistic Schrodinger equations in the presence of the QCD potential. These functions 
are defined by 



\- VgcDir 

nit 



mt 



+ VQCD{r) - {E 



G{^,E) = 
F\^,E) 



6'i 



k r3/ 



and 



G{p,E) 
p'F{p,E) 



d^xe-'^''G{x,E), 



(2.7) 
(2.8) 

(2.9) 
(2.10) 



One may obtain the Green's functions either by first solving the Schrodinger equations in 
coordinate space and taking the Fourier transforms of the solutions @], or by solving the 
Schrodinger equations directly in momentum space |Q. 

Various known 0{as) corrections to threshold cross sections can also be expressed in 
terms of the above Green's functions. In Ref. ||T^ the following functions have been defined, 
which are solely determined from QCD, to represent various independent corrections. We will 
pursue the same conventions in this paper. 



^ip,E) 



:i-C^a,(mi)/7r) p F\p,E) 



(1 - 2GFasimt)/TT) mt G\p, E) ' 
Rey?, = Im^j, 



(2.11) 
(2.12) 
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d3q 



- CF-47ra^(/iB) Pr. 



1 



|q-pt| 



2Im 



G{q,E) 



1 



Pr(q-Pt 



'2 



(27r)3 



|q-Pir iPtI |q-ptl 



2 Re 



G{q,E) 



G{p, E) 



(2.13) 



(2.14) 



Here (7^ = 4/3 is a color factor; p = |pf| and g = |q|; /x^ is taken to be 15 GeV in our analyses. 
Eqs. ( |2.13|) and ( p.l4| ) differ slightly in their integrands from those defined in Ref. [l^.Q The 
differences, however, can be regarded as higher order corrections which are beyond the scope 
of our analysis. 



3 Qualitative Description of Final-State Interaction Ef- 
fects 

In this section a qualitative understanding of the effects of final-state interactions is explained. 
It is based on the classical picture that t and h (t and h) attract each other due to their colour 
charges. We will see that all the following qualitative features match well with our numerical 
studies presented in section 5. Moreover, the following argument will help interpreting the 
formulas derived in the next section. 



3. a Top Momentum Distribution 

Let us first consider the effect of final-state interactions on the top-quark momentum (|pt|) 
distribution, where pt is reconstructed from the hW^ momenta at time r — > cxd. The energy of 
a ti system before its decay is given by 

8 = Et + E-t + V{\rt-Vt\), (3.1) 

where 

Et = mt + ^, E, = mt + ^, V{r) = -Cp^. (3.2) 

Suppose t decays first and let a ~ {asmt)~^ ~ (m^rf)"^/^ be the typical distance between t and 
t at the time of t decay. Then, just before the decay, the momenta of t and t are given by 



\Pt\ = \Pt\ ~ \/mt [S - 2mt + \V{a)\]. (3.3) 

Their order of magnitude is asfrit ~ (Bohr radius)"^. If it were not for the final-state inter- 
actions between i and b, and if i, b and traveled as free particles, the above momentum 
would be transferred to the bW~^ system at time r — > 00: 



Pb + Pw+\t^oo Jmt[£ - 2mt + \ V{a)\]. (3.4) 



* In Ref. IT^] the Coulomb propagator between h and i (or between h and t) together with the color charges 
are replaced by the QCD potential between two heavy quarks by hand, — CF-47ras/|q— ptp Vqcd{\(1— Pt|)- 
One advantage of Eqs. (2.13) and (p. 14 ) is that one may convert them into one-parameter integral forms by 
explicitly integrating over dilq |jl^, |13|. 
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Figure 1: Attractive Coulomb force between t and b (from t decay). The momentum transfer 
6ph = —Spt due to the attraction is indicated by thick arrows. 



Taking into account the final-state interaction (Coulomb interaction) between t and b, the 
energy of the ibW^ system is given 



£ = Eh + Ew++Et + Vi\n-rt\), (3.5) 



Et, = |p,|, Ew+ = Vp'h^+ + M^. (3.6) 

As depicted in Fig. |I], the tracks of b and t are deflected due to the attraction between the two 
particles, which lose kinetic energy as b flies off to infinity at the speed of light. The classical 
equation of motion is given by 

dpb dpt d 

-1- = = -i^^(|r6-rf|). (3.7) 
ar ar orb 

Substituting the free particle solution rj, = vr+ro on the right-hand side and noting |rb— rj | ~ r, 
we may estimate the size of the momentum transfer due to the attractive force as[| 

\5pb\ = \Spt\ ~ \Virrmn)l (3.8) 

where the minimum distance between b and i is denoted as rmin- Typically rmin ~ ct, and we 
find from Eq. (|3.8|) that the effect of final-state interactions on the \pt\ distribution is of the 
order of \6pt\ ~ a^nit. Obviously the effect is to reduce |pt|. 

3.b Forward-Backward Asymmetric Distribution 

Next we consider the cos 9te distribution of the top quark. {9te denotes the angle between t and 
e~ in the ti cm. frame.) It has been known that a forward-backward asymmetric distribution 
of the top quark is generated by the final-state interactions [l^, |l^ . We describe its mechanism 
here. 



* Here we neglect the interaction of i (b) and the magnetic field generated by b (t). This approximation is 
justified in the case of our interest; see section 4. 

t It corresponds to solving the equation of motion by a series expansion of a^. 
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Wj-^^< b 

Figure 2: Typical configurations in the decay of i witli definite spin orientation. Transverse 
W~ (Wrf) tend to be emitted in the direction of the t spin orientation, while longitudinal W~ 
(W^) are emitted in the opposite direction due to helicity conservation. For mt ~ 175 GeV, i 
decays mainly to W£, hence b is emitted more in the t spin direction. 



We consider the case where i decays first and examine the interaction between t and b. 
The t and i pair-produced near threshold in e^e~ collisions have their spins approximately 
parallel or anti-parallel to the e~ beam direction and the spins are always oriented parallel to 
each other. In fact in leading order the polarization vector of the top quark is given by 



P = C°(x) n||- On the other hand, the decay of i occurs via a V—A coupling, and b is emitted 
preferably in the spin direction of the parent t, see Fig. More precisely, the excess of the 
6's emitted in the t spin direction over those emitted in the opposite direction is given by k 
defined in Eq. ( |2.6| ). Now suppose t and t have their spins in the ny direction. Then b will be 
emitted dominantly in the n|| direction. One can see from Fig. |^(a) that in this case t is always 
attracted to the forward direction due to the attractive force between t and b. The direction of 
the attractive force will be opposite if t and i have their spins in the — ny direction (Fig. 0(b)). 
Thus, polarized top quarks will be pulled in a definite (forward or backward) direction, and we 
may expect that a forward-backward asymmetric distribution of the top quark ~ k,C^^{x) cosOte 
is generated by the final-state interaction. 

Incidentally, a forward-backward asymmetric distribution of the top quarks is also gener- 
ated by the interference between the S'-wave and P-wave ti pair-production amplitudes It 
is formally a quantity of the same order as the final-state interaction near threshold, since it 
arises as an 0{(3) = 0{as) correction to the leading spherically symmetric distribution. Inter- 
estingly, we find that each cosOte distribution has quite a different physical explanation for its 
generation mechanism, although it is a common feature that both originate from the interplay 
between QCD and electroweak interactions. 

3.C Top Quark Polarization Vector 

The Coulomb attraction between t and b also modifies the top-quark polarization vector |]TU 



As we have seen previously, in the i decay, b tends to be emitted in the direction of the parent 
?s spin direction (Fig. We then find from Fig. ^ that if the t and i spins are oriented in 
the nil direction, t will be attracted to the forward direction due to the attraction by 6, and 
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t 



t 




1^1 



A 



1+ 



b 



W 



(a) 



(b) 



Figure 3: Attractive force between t and h when the t and i spins are oriented in the (a) 
ny direction, and in the (b) — ny direction. The momentum transfer 5p;, = — 5pf due to the 
attraction is indicated by thick black arrows. 



oppositely attracted to the backward direction if the t and t spins are in the — n|j direction. 
This means that in the forward region (cos 6'ie — 1) the number of t's with spin in ny direction 
increases whereas in the backward region the number of those with spin in the opposite direction 
increases. Or equivalently, the ny-component of the top-quark polarization vector increases in 
the forward region and decreases in the backward region. We may thus conjecture that the 
top-quark polarization vector is modified as ST* ~ n cos 6te ny due to the interaction between t 
and h. 

3.d Z+ Energy-Angular Distribution 

Finally let us examine the effect of the Coulomb attraction between h and i on the energy- 
angular distribution in the semi-leptonic decay of t. The 6-quark from t decay will be attracted 
in the direction of i due to the Coulomb interaction between these two particles. We show 
schematically typical configurations of the particles in the top-quark semi-leptonic decay in 
Fig. ^ It can be seen that if the probability for h being emitted in the i direction increases, 
correspondingly the probability for particular energy-angular configurations increases. These 
configurations are either "i?; is small and emitted in — direction" or "i?; is large and 
emitted in p^ direction". 
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t t 




l*^ ' ^ V "Ej small, /+in -p^ direction 

V ^ '■ ^ /"^ "Ej large, /"^ in direction" 

Figure 4: Typical configurations of the particles in semi-leptonic decay of t when the 6-quark 
is emitted in the i direction. Due to the boost by W~^, the energy-angle correlation of l'^ will 
be either "i?/ is small and emitted in —pt direction" or "E'i is large and emitted in pt 
direction". 




Figure 5: Born-type diagram for the process e+e tt ^ hl^uhW . 

4 Z+ Energy- Angular Distribution 

We present the formulas for the charged lepton energy-angular distribution in the decay of a 
top quark that is produced via e^e~ —>■ tt near threshold. In the following, xi = 2Ei/mt and fli, 
respectively, denote the normalized energy and the solid angle of the charged lepton as defined 
in the rest frame of the parent top quark. For simplicity we neglect the decay of W~ in our 
calculations. 

4. a Factorizable Part 

It is well-known that the contribution of the Born-type (reducible) diagram (Fig. ^ to the 
differential distribution of t and has a form where the production and decay processes of the 
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top quark are factorized: 



dcTBornje^e ^tt ^ hl^vhW ) _ daBomje^e tt) ^ J_ dTt-^U+u CP Bom) . X 

d^Ptdxidili d^pt dxidVli 

The above form holds true even including 0{as) = 0{(3) corrections to each vertex and prop- 
agator in the Born-type diagram. Here, Tt ||2^, and dT t^bi+^iT* Bom) /d^idfli p^, |16 



re- 



spectively, are the width of a free top quark and the charged lepton energy- angular dis- 
tribution in the decay of a free polarized top quark, both including the corresponding full 
0{as) corrections. Near threshold, the top-quark production cross section for longitudinally 
polarized e^e~ beams is given by |T^, |T9[ 



d(7Bomie^e ^ tt) da% ,^ , ^ , n i 

d^t ^d^t"" [1 + 2CFB(x)</'K(p,^)cos^te] , (4.2) 
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d^pt^ Anmf a-^e+Pe-)(ai + X«2) (^1 ' j G{p,E) . (4.3) 

Here, a is the fine structure constant, and Nc = 3. T-'Bom = "Pw^W +'P±n± + 'PNnN represents the 
polarization vector of a top quark produced via the Born-type diagram (Fig. ^ near threshold. 
The components are given by 



V\\{pt,E,x) = Cf(x) + q}(x)<^R,(p,i^) cos^te 
'P±{pt,E,x) = C±{x) Vb^P^E) sin Ote 
^N(Pi,^,x) = Cn(x)v?i(p,^) sin 6*46. 



(4.4) 
(4.5) 
(4.6) 



Let us review briefly how to derive the factorized form of the differential distribution, 
Eq. ([4.1|) . First, in calculating the fully-differential cross section da/d^r,{bl'^i'bW~), one may 
replace the top-quark momentum by an on-shell four- vector as 



P^={p't,Pt) - pf = Wp't+mlpt) (4.7) 

in vertices and propagator numerators (but not in propagator denominators). The replacement 
is justified because near threshold relevant kinematical configurations are determined by 

p°-mt^almt, |pt| ~ a^m* (4.8) 

so that the replacement induces differences only at (9(0^), and also because we will not be 
concerned with the p°-dependence of the cross section. Then, one may use the following iden- 
tity to factorize the spinor traces that appear in the fully-differential cross section into their 
production and decay parts: For an arbitrary 4x4 spinor matrix G and for a four-vector 
satisfying pi = m^, 

^t + m ^ ^t + m ^ :^t + m i- 'pj^ ^ 

2mt 2mt 2mt 2 ' ^ ' 
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where the four-vector and the constant C are determined from G and via the relation 

2 

provided and satisfy s ■ pt = V ■ pt = 0. 
One may also factorize the phase-space as 

74 

(i$5(7* ^ bl+ubW) = d^3{t* btv) d^2{t* bW-), (4.11) 

(27r)4 

d^s{t* ^ = (^i^) dxmidplrd(t)u, (4.12) 

where is the invariant-mass-squared of and 0^ denotes the azimuthal angle of b around 
/+ in the top-quark rest frame. Then one integrates over dp^, d(j)bh d^2(bW~) and dp^/{27i); 
the integration over is trivial since we use the narrow-width approximation for W~^] the 
integration over d(j)i,i is also trivial since the fully-differential cross section is independent of 0^; 
the integration over the bW~ phase-space merely replaces the bW~ wave functions by r^; the 
p° / (27r)-integration is straightforward. 



^t + ^t ^ ^t + m 1- 



2mt 



2mt 



(4.10) 



4.b Final- State Interaction Corrections 

Corrections due to the final-state interactions (rescattering corrections) that originate from 
the irreducible diagrams (a)-(d) in Fig. ^ are important particularly in the threshold region. 
In fact their contributions are counted as 0{as) = 0{P) corrections to the leading threshold 
enhancement [|10|. We calculate the effect of each diagram on the energy-angular distribution. 
We chose Coulomb-gauge for the QCD part in our calculations. 

The contribution of diagram (a) (exchange of one Coulomb-gluon between t and b) can be 
regarded, after integrating over the bW~ phase-space, as a correction to the production process 
of t. Thus, the production cross section (i(TBom(e"'"e~ — > ti)/d^pt and the polarization vector 
'^Born of the top quark receive corrections by this diagram, whereas the decay distribution 
dTt^bi+^iT'Born) / dxidfli remains unaffected (except for the modification of 7-*Born)- In fact the 
contribution of this diagram can be incorporated by the following substitutions in Eq. ( [4.11 ): 



C?0-Born(e+e" 



tt) 



da 



Bornle^e 



tt) 



d^Pt 



d^Pt 



X (1 + 5a) 



Born 



VBorn + SVa (4.13) 



with 



Sa 



Ipl {p, E) + KC?\ij^ {p, E) cos Ote 



I -{Clf\Ki,^{p,E) COS ete-ii\\. 



(4.14) 
(4.15) 



The derivation of the formula goes as follows. Since the relevant kinematical configuration 
lies in the soft gluon region, we can use soft-gluon approximation and factor out the part that 
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(c) (d) 

Figure 6: Diagrams for final-state interactions for e'^e~ — > — > bl'^ubW~: (a) Coulomb-gluon 
exchange between t and b, (b) Coulomb-gluon exchange between i and b, (c) Coulomb-gluon 
exchange between b and b, and (d) transverse-gluon exchange between b and b. 



11 



depends on the loop-momentum (the propagators of the gluon, b, t and t together with 
the loop-integral) outside the spinor structure, while the remaining part is similar to the fully- 
differential cross section of the Born-type diagram. The latter part is factorized as before. Due 
to the soft-gluon factor (the factor pulled outside), at this stage we may interpret that the 
production cross section and polarization vector of top quark get corrections that depend on t 
and h momenta. Integrations over and c/^f,; are the same as for the Born-type diagram. 
For integrations over d(^2{pW~), dp^/{27T) and dq^/{2TT), we follow the method described in 
Ref.[|Tg, Appendix D. We are thus led to Eqs. (^:T3D -(p^. 

Diagram (b) (exchange of one Coulomb-gluon between i and b) in Fig. |^ gives a correction 
that connects the production and decay processes of the top quark. In fact one may incorporate 
the contribution of this diagram by multiplying Eq. ( [4.1|) by a factor [1 + ^{p, E,xi, cos 9^)], 
where 



^{p, E, xi, cos 6it) = Cf ■ AirasifJ^B) 



G*iq,E) 

|q-ptP '^^ [G*{p,E) 



Re 



2n 



2-K 



|q-ptl 



n6-(q-Pt) + «e 



(4.16) 



fi;, denotes the unit vector in the direction of h. After integration over dQcid4'bh oiie niay reduce 
the expression to a one-parameter integral form as 



^{p, E, xi, COS 9it) = Cf ■ 4:Tcas{nB) / dqlw^Re ^ 

JO I LG* 



G *{q,E) 



+ Im 



G*{q,E) 



G*{p,E) 



(4.17) 



with 



W- 



47^2 p2 _ q2 

1 q 



I 4vr2 \p2 _ q2 



^ {^(4- l)cosh-^ |z+| -0(£- l)cosh-^ (4.18) 

■ [tvO ((1 + z+){l + z-)) - 6(1 - 4) cos-^(;z+) - ^(1 - z'i) cos"^(z_)} , 

(4.19) 



z± = 



q p ± cos 9wi cos 9it 1 + y 2y ^ Pt • P« 

= ^ ' cos6'h-/ = - -, cos6'it = -^ — — 

sm9wism9it l-y xi{l - y) \Pt\\Pi 



\Pt\\Pi\ 



(4.20) 



Here, 9wi represents the angle between and /"*" in the t rest frame (given as a function of 
xi); 9it represents the angle between t and in the tt cm. frame[| (0 < 9wi,9it < tt). The 
inverse functions cosh~^ and cos~^ in the above formulas take their values within [0, oo) and 
[0,7r], respectively. 9{x) is the unit step function. It is understood that the principal value 
should be taken in the integration of the w^-term as p — > g. 

We derived Eqs. ( [4.16| ) and ( [4.17| ) in the following manner. As in diagram (a), we used the 
soft-gluon approximation and factored out a soft-gluon factor (loop-integral of the propagators 
of gluon, b, t and t). The remaining part is same as the fully-differential cross section of 
the Born-type diagram except for the t and i propagators and Green's functions, which is 



* Within our approximation, there is no distinction between cos On and — cos 6it, where 9it denotes the angle 
between i and in the t rest frame. 
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again factorized. This time, however, the correction cannot be interpreted as associated with 
the top-quark production process since the soft-gluon factor depends on the b momentum. 
The integrations over dp'^y and d^2{bW~) are the same as those for the Born-type diagram. 
The function ^ given in Eg. ( |4.16D is essentially the soft-gluon factor integrated over (ip°/(27r), 
d'^q/(27i)^ and d(f)bi. Finally, to derive Eq.( [4.17| ) from ( |4.1(j| ), it is simpler to integrate over dfl^ 
before ^0^. 

Two noteworthy properties of ^{p, E,xi, cos On) are: (1) its /"""-angular dependence enters 
only through cosdu and is independent of the angle from the e~ beam direction or from the 
top-quark polarization vector, and (2) it is purely determined by the QCD interaction and free 
of the coupling parameters of electroweak interactions (except y). 

As a non-trivial cross-check of the formula (|4.17|) , we integrated ^ x T^^dTt^hi+u/dxid^li 
over the lepton energy- angular variables / dxidQi analytically and reproduced the one-parameter 
integral formula [10, 13| for the final-state interaction correction (from diagram (b)) to the top- 
quark three-momentum distribution. 

The contribution of diagram (c) (exchange of one Coulomb-gluon between b and b) van- 
ishes within our approximation. We show it in steps. Using soft-gluon approximation, the 
contribution of this diagram to the cross section can be written as 

xzCp-inasJ (0 [D{pt) + D{pt)] [D*{pt + k) + D*{p,-k)] 
xGip,E)G*i\pt+k\,E) 



A;*^ — fif, ■ k — ie — + fij-k — ie |kp 
+complex conj., (4-21) 

where k^^ is the gluon momentum and 

= rf-.n.,-p;/2.„. + .r,/2 (^■^2) 

denotes the non-relativistic top-quark propagator. T is the contraction of hadronic and leptonic 
tensors resulting from the spinor traces after a soft-gluon factor (the second and third lines) is 
taken out. It coincides with the fully-differential cross section of the Born-type diagram except 
for the t and i propagators and Green's functions; hence T is real. Next we integrate over 
(ip°/(27r), dk°/{2n) and d^k. Let us define 

/(k, hb, hi) 

(4.23) 



Noting that the contribution of diagram (c) to the cross section comes solely from the gluon 
momentum region 

k\ |k| ~ a^^mt, (4.24) 
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(e) 



(f) 



Figure 7: Diagrams for real-gluon emission process e~^e tt ^ hl^uhW g. 



we keep |k| within this region. Then we may substitute G*{\pt + k|,£') —>■ G*{p,E) in Eq. 
( [4 .211) since the difference is higher order, and find 



It is easy to see that / dfly^ I(k, hb, hi) is real using the symmetry of dfl^ under k — >• — k. Thus, 
we conclude dac/d^ptdxidfli = in our approximation. 

In fact the same proof can be applied to show quite generally that the contribution of 
diagram (c) vanishes at 0{as) provided one calculates a cross section where the top-quark 
energy is integrated out; for example, the top-quark three-momentum distribution. This is 
no longer the case when one considers a cross section that depends explicitly on the top-quark 
energy; for example, the top-quark four-momentum distribution. Then the diagram in question 
does contribute. 

One can also show in a similar way that the kinematical regions Eq. ( [4.24|) in diagrams 
(a) and (b) do not contribute to the cross section da / d^ptdxidQi and that only the gluon 
momentum region where |k| ~ a^mj ^ Ft is relevant. We took advantage of this fact in 
deriving Eqs. (^13D-(^l7D. 

It can be shown using similar techniques that the contribution of diagram (d) (exchange of 
one transverse gluon between b and b) gets canceled when it is added to that of the corresponding 
real-gluon emission diagram (interference of diagrams (e) and (f) in Fig. 0). This cancellation 
is consistent with the same cancellation that was found in the calculation of the top-quark 
momentum distribution [^. The contribution from each of these diagrams comes from the 
gluon momentum region |k| < a^m^ and is in fact logarithmically divergent due to a soft-gluon 
singularity. 




+ complex conj. 



(4.25) 
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The cancellations of the final-state interaction corrections at the various levels of inclusive 
cross sections are summarized in section 7. 



Now let us compare our formulas and the argument given in the previous section. In 
subsections 3.b and 3.c, the parity-violating nature of the electroweak interactions in top pro- 
duction and decay played an essential role, while this was not the case in subsections 3. a and 
3.d. We find that correspondingly the cos6'fe term of Sa and 67^a contain electroweak coupling 
parameters (through C||), while the symmetric term of 6a and ^ are independent of these elec- 
troweak parameters. More precisely, the cos 6te term of Sa and S'Pa have the forms anticipated 
in subsections 3.b and 3.c, respectively, if the function ip^{p,E) is positive. Indeed, the nu- 
merical evaluation in Ref. [0 shows that tp^ {p, E) ^ holds in the entire threshold region. 



in d'Pa can be understood within our previous 



Besides, an additional coefficient 
argument in the extreme cases C^^ = ±1. Namely, if the top quark is 100% polarized, there 
will be no contamination from the opposite spin so that the correction should disappear. It 
may be interesting to note that the final-state interaction corrections to the polarization vector 
vanishes for the ideally polarized top quarks, Ci? = ±1. 



4.C Formula Including Full 0{as) Corrections 

In summary, the energy-angular distribution of including full 0{as 
be cast into a form 



0{P) corrections can 



da{e+e- ^ tt ^ bl+ubW') da{e^e 



with 



d^Ptdxidfli 



da{e'^e tt) 



d^Pt 



tt) 1 dVt-,u+.{V) n , c\ 

— X ^ — X (1 + 

Vt dxidVti 



d^Pt d^pt 



(4.26) 

(4.27) 
(4.28) 



The above distribution is obtained as the sum of the cross sections for e"'"e~ — >■ — > hl^vhW~ 
and e^e~ tt bl'^ubW' g. An independent emission of a gluon from the t or t side has been 
included in Eq. ( [4.1|) , while the interference of both has been incorporated in conjunction with 
the final-state interaction diagram (d) in subsection 4.b. 

Experimentally the top-quark four-momentum p^* will necessarily be reconstructed from 
the bW~ system in the study of the distribution. For the case with a gluon in the final state, 
we assign the "top-quark momentum" as 



Case (A) : pt = {pe- + Pe+) - {Pl + PW- + Pg) if {Ph + PW' + Pg 



l2 _ ^2 < 



Case (B) : pt = {pe~ + Pe+) - (pi + Pw 



if {pi + pw- 



t ~ 



t ~ 
mt^t 



mtV 



(See, however, the discussion in section 7.) In other kinematical configurations the cross section 
is suppressed. Experimentally there will be a corresponding cut in the hW~ invariant mass. 
If both conditions in (A) and (B) are satisfied simultaneously, the gluon should necessarily 
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Top Momentum Distribution [pb/GeV] 
0.05 I 
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Figure 8: Top-quark momentum distributions da/d\'pt\ obtained from Eq. ( [4.13| ) for various 
cm. energies measured from the lowest-lying resonance, /S.E = y/s — Mis- Solid lines and 
dashed lines represent the distributions with and without the final-state interaction between t 
and h {5a), respectively. We set Pg- = 0. The left figure corresponds to the coordinate-space 
calculation, the right one to the momentum-space calculation. 



be soft, and there will be no difference within our approximation between the cross sections 
corresponding to the above two assignments of the top-quark momentum. 

One comment is in order here. In defining the "rest frame" of a top-quark in the study of 
distribution, one may use either or (defined in Eq.( [4.7|) ). The difference of the cross 
sections based on the two definitions is of 0{a1) which is beyond the scope of our approxima- 
tions. Thus, the cross sections defined in both definitions should be measured in experiment 
and compared. It will serve as a cross-check for the stability of our prediction. 



5 Numerical Results 

In this section we examine the effects of the final-state interaction corrections on the energy- 
angular distribution numerically, and compare the results with the qualitative argument given 
in section 3. The numerical results are obtained using both the coordinate-space approach 
developed in Refs. p|, |[ and the momentum-space approach developed in Refs. 0, Con- 
ventionally these two approaches have been used independently by different groups, and this is 
the first time to make a direct comparison of the cross sections calculated in both approaches. 
Some of the produced results are slightly different. We set rrit = 175 GeV, as{Mz) = 0.118, 
Pe+ = and a = 1/128 in all our analyses. 

We first examine the contribution of diagram (a) in Figs. ^ (Coulomb interaction between t 
and b) as given in Eqs. ( ^.13| )-( p:.15D . Shown in Fig. ^are the top-quark momentum distribution 
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dij/d\'pt\ for various cm. energies measured from the fowest-lying resonance, Q A = ^/s — Mls. 
This is calculated from the top-quark production cross section Eq. ([4.13|) . As expected, the 
top-quark momentum is reduced. The effects are half in magnitude as compared to the final- 



state interaction corrections given in Refs. ||10|, since only the interaction between t and h 
is included here. 

We show the angular distribution of the top-quarks in Figs. |^. It can be seen that the final- 
state interaction increases the top-quark distribution in the forward direction for P^- = 1. This 
is consistent with our argument in subsection 3.b since in leading-order approximation t and t 
have their spins aligned perfectly in the ny direction for this e~ polarization. Oppositely we see 
that the final-state interaction decreases the top-quark distribution in the forward direction for 
Pg- = — 1. We note that the top quark has a natural polarization ~ — 0.4n|| for unpolarized 
e^e~ beams. Hence, the sign of the correction is the same as in the P^- = —1 case. Also 



we show corrections to the top-quark polarization vector in Fig. |Ty. Although the qualitative 
behaviour meets our expectation, the magnitude of the correction is rather small. Note that 
(57-*a vanishes for Pg- = ±1 since C^^ = ±1. 

Next we investigate the irreducible (non-factorizable) correction that stems from the 
Coulomb interaction between i and b. The correction factor ^{p, E,xi, cos On) given in Eqs. 



( [4.17|) -( ^.20 ) depends on four parameters, two of which specify the lepton configuration: xi 



and cos 6it. Therefore, we will examine the dependence of ^ on these two parameters for several 
(p, AE) combinations. We fix the top-quark momentum p to be the peak momentum of the dis- 
tribution for each AE, hence its values are slightly different for the two numerical approaches; 
see the top-quark momentum distributions in Fig. |[ Shown in Fig. |Tl| are 3-dimensional plots 
of ^ as a function of xi and cosOu- One can see that in the figures ^ takes comparatively large 
positive values for either "small xi and cosO^ ~ —1" or "large xi and cosOn ~ Oppo- 
sitely, in the other two corners of the x^-cos On plane ^ becomes small or becomes negative for 
smaller AE. The typical magnitude of ^ is 10-20%, which would be a reasonable size for an 
0{as) = 0{/3) correction. This behaviour holds also true for ,^ at p off the peak of the top- 
quark momentum distribution. These features of the correction factor C, are consistent with our 
qualitative argument in section 3. 

We made a cross-check of our numerical results for C, by numerically integrating ^ x 
T^^dTt^i^i+i^/dxidfli over the lepton energy-angular variables / dxidQi and comparing to the 



final-state interaction corrections to top-quark momentum distribution given in Refs. ||I0| , |19 

It is seen that in Figs. || and ^ the coordinate-space approach and the momentum-space 
approach produce slightly different results. In particular the normalization of the cross sections 
differ at lower cm. energies, ^ AE < 2 GeV, whereas the differences decrease at higher 
energies. The cause of the differences can be traced back to the different short-distance QCD 
potentials employed in the two approaches. As we will discuss in section 7, this difference is 
formally counted as higher order beyond our approximation, and at present it should be taken 
as an uncertainty of the theoretical prediction. 



* Mis is defined as the real part of the position of the lowest-lying resonance pole in the complex energy 
plane. The energy measured from Mis is more convenient than the energy measured from the threshold 
{E — ^/s — 2mt) when we compare different potentials in the literature. 
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a) Coordinate-space approach results 
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b) Momentum-space approach results 



Figure 9: Top-quark angular distribution da / dcosdte obtained from Eq. ( |4.13| ) for various 
electron polarizations P^-- Solid lines and dashed lines represent the distributions with and 
without the final-state interaction between t and b (5a), respectively. 
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a) Coordinate-space approach results 
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b) Momentum-space approach results 



Figure 10: Final-state interaction correction to the ny component of the top-quark polarization 
vector, SVa, for P^- — -1-0.8 (sohd line), (dash hne), —0.8 (dot-dash line). 
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Figure 11: Three-dimensional plots of ^ as a function of xi (x-axis) and cosdte (y-axis) at the 
peak momenta of the |pf| -distribution in Fig. ^ AE = GeV and AE = 5 GeV. 
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6 Observable Proper to Top Decay Processes 



As we have seen in the previous section, the final-state interactions affect the distribution 
in top-quark decays. The correction factor ^ depends on the kinematical variables of both the 
top quark and /''', and destroys the factorization of the cross section Eq. ( [4.1| ). In this section 
we define an observable which depends only on the top-quark decay process {dTf-^bi+u/dxidfli 
of a free polarized top quark) It is a differential quantity dependent on the energy- angular 
variables. 

From Eq. ( 4.17 ) one sees that the correction factor ^ is invariant under the simultaneous 
transformations of the angular variables 



cos 6 



Wl 



cos 6 



Wh 



COS 9it - cos 9it, 



(6.1) 



since z± is invariant. This invariance may be understood as follows. The soft-gluon factor, 
which represents the final-state interaction in diagram (b) (Fig. |^), does not depend on the 
momentum as long as the momentum is kept fixed. One sees that accordingly the integrand 
of Eq. ([4.1 6|) is independent of the energy and angle. The dependence on these variables 
enters only through the phase-space integration over (pu for a given configuration. There- 
fore, one may reverse the momentum in the t rest frame without affecting the phase-space 
integration, thereby keeping the whole function ^ also unchanged. The above transformations 
Eq. ( |6.1|) are essentially this reversal of the momentum. (Due to the form of the integrand, 
there are extra degrees of freedom for the transformation of the direction, see below.) 

Using Eq. ( [4.20|) , the above transformation can be written as a transformation of the 
lepton energy and angle as 



l + y_ 

, y 



n. 



(6.2) 



Here, hi = Pi/\pi\ denotes the unit vector in the direction of in the top-quark rest frame. 
The choice of fi; for flipping the sign of cosO^ {cosOu — > — cos 9it) is not unique. We represent 
by fi; an arbitrary one of those choices. The production cross section of the top quark is not 
affected by this transformation since the top-quark kinematical variables are not involved. The 
important point is that neither the flnal-state interaction correction is affected by it. 

Now, using this invariance, we flrst construct a quantity that has the simplest structure 
from the theoretical point of view, and afterwards we present an improved quantity that will 
be more useful for practical purposes. Let us deflne 



A{xi,hrV,h'i-V) = 



da{e+e-^tt bl+ubW- 
d^PtdxidQi 



'da{e+e-^tt^bl+ubW- 
d^PtdxidQi 



(6.3) 



* This property is true only up to 0{as) — 0{/3) corrections and may be violated by yet uncalculated 0{aj.) 
corrections. 
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The production cross section and the correction factor ^ cancel in the numerator and denomina- 
tor. As a result, this quantity is independent of the top-quark momentum and is determined 
only from the free polarized top-quark decay cross section. In fact, using Eq. ( 4.26 ), we find 
that 



A{xi,hrV,h[-V) 



dT 



t^bl+u[ 



dxidVti 



dT 



dxidQi 



(6.4) 



holds up to (and including) 0{as) = 0{P) corrections. Note that the polarization vector 
7-*, which specifies the decay distribution, includes the correction induced by the final-state 
interaction between t and b; see Eq. ( [4.2 



Since we take a ratio of differential cross sections in the definition of A in Eq. (|6.3| ), this 
quantity would suffer from a large statistical error experimentally. Meanwhile, this quantity is 
predicted to be dependent only on a few variables theoretically. (See Eq. (|6.4] ).) It means that 
we may integrate out the irrelevant kinematical variables before taking the ratio and reduce its 
statistical uncertainty.^ For instance, we may choose nj = — ri/ and define 



A{xi, a) = 



1 d^ptdQi 6{hrV-a) 


' da{e+e- tt ^ hl+vhW-) 




d^PtdxidQi 


J d^ptd^i 6{hrV + a) 


'da{e+e- ^ tt ^ bl+vbW-) 




d^PtdxidQi 



(6.5) 



Here, the top-quark polarization vector in the delta functions should be evaluated as a func- 
tion of pt according to Eqs. (|4.4| )- (|4l^ ), ( |4.15| ) and ( |4.28| ). The numerator and denominator, 
respectively, depend on two external kinematical variables and all other variables are integrated 
out. 

Again using Eq. (^4.26|) , one finds that A is determined solely from the free polarized top 
decays: 



A{xi,a) 



dV 



dxidVti 



h,7>=a • 



dT 



dxidVti 



Xl^Xj 

xilT>=~a 



(6.6) 



t Consider a ratio of certain physical quantities which depend on a set of kinematical variables 

K 



where 4>i denotes a point in the phase space, and (p^ is obtained by a transformation from it, — 
Whenever the ratio R takes the same value in a subspace U of the whole phase space, we may take sum over 
the subspace before taking the ratio: 



R 
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This is a general formula that is valid even if the decay vertices of the top quark deviate from 
the standard-model forms.| We see that the quantity A{xi, a) preserves most of the differential 
information^ contained in dTt^hi+^/dxidQi. 

7 Discussion 

In this section we discuss three different issues relevant to our work. These are: the difference 
between the coordinate-space and momentum-space potentials, the mis- assignment of the top- 
quark momentum, and the disappearance of the final-state interaction corrections at the various 
levels of inclusive cross sections. 

As we have seen in section 5, our numerical results obtained from the coordinate-space 
calculations and those obtained from the momentum-space calculations differ slightly, although 
all the qualitative features are common. The difference can be traced back to the difference in 
the short-distance part of the QCD potentials used in the two approaches. 

Let us remind the reader how each potential is constructed (in the short-distance regime). 
The large-momentum part of the momentum-space potential Vjkt 0? i] is determined as fol- 
lows. First the potential has been calculated up to the next-to-leading order in a fixed-order 
calculation. The potential is then improved using the two-loop renormalization group equa- 
tion in momentum space. On the other hand, the short- distance part of the coordinate-space 
potential Vsfhmn |@] is calculated by taking the Fourier transform of the fixed-order potential 
in momentum space, and then the potential is improved using the two-loop renormalization 
group equation in coordinate space. Thus, the two potentials are not the Fourier transforms 
of each other. Only the leading and next-to-leading logarithmic terms of the series expan- 
sion in a fixed MS-coupling are the same for the two potentials. The difference begins at the 
next-to-next-to-leading order terms. (The non-logarithmic term in the two-loop fixed-order 
correction.) 

To make a clear comparison, the two potentials are Fourier transformed numerically and 
we examine their difference both in coordinate space and in momentum space. We show the 
effective charges defined as 



.(1/r) = (-Cp/ry'Vir), a^^Q) = (-47rC^/g') V^(g) (7.1) 



in Fig. [T^, which clearly demonstrates that there is a non-negligible difference between the 
potentials. The oscillatory behaviour of Vsfhmn in Fig. |l^(b) is an artifact due to the dis- 

Note that, quite generally, energy-angular distributions of from free polarized top quarks have a form 

— — = Fo{xi) + {nrr)Fi{xi). 

dxidili 

§ According to its construction, it satisfies a relation 

A{xi,a) ■A{x'i,-a) = 1. 
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Figure 12: The effective charges (a) acoor(l/''^) and (b) amom(p) defined via the potentials in 
coordinate space and momentum space, respectively. Vjkt [li and Vsfhmn represent the 
potentials used in the momentum-space approach and the coordinate-space approach, respec- 
tively. 



continuity in the second derivative of the coordinate-space potential, which is located at the 
continuation point of the perturbative potential to a long-distance potential. As already stated, 
the coordinate-space potential Vsfhmn follows the form required by the two-loop renormaliza- 
tion group equation in the short-distance region, whereas the momentum-space potential Vjkt 
follows the form required by the two-loop renormalization group equation in the large momen- 
tum region. 

In principle we may reduce the difference by including the two-loop finite correction and 
invoking the three-loop renormalization group improvement [BDl. We shall not do so in this 



paper, because there are a number of other corrections of the same order of magnitude which 
are not calculated yet. We will study the difference of the two approaches in more detail in a 
forthcoming paper. 

In subsection 4.c, we assumed a perfect assignment of the top-quark momentum in cases 
(A) and (B) for defining the distribution formula, Eq. ( |4.26| ). In real experiments, however, 
a mis-assignment of the top-quark momentum will be inevitable whenever there is real gluon 
radiation in the final state because of the typical jet-clustering algorithm that will be used. For 
instance, when a gluon is indistinguishable from a 6-jet, the top-quark momentum reconstructed 
by clustering may be off-shell; rather grouping the gluon on the other side (with b) would result 
in an on-shell momentum.^ Ref. studied how this mis- assignment alters the top-quark three- 



momentum distribution near threshold and found that the correction is less than a few percent; 
the clustering algorithm assumed in that paper, however, is somewhat unrealistic. The effect 



* There is no ambiguity in assigning the gluon to the production of tt since real gluon radiation in the 
top-quark production process is suppressed near threshold. 
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T 1 

mclusiveness 


Diagram (a)/(bj 


Diagram (cj 


Diagram (d) + (e)*(f) 


fully differential 


+ 


+ 


+ 


top energy integrated 


+ 


vanish 


vanish 


O-tot(s) 


vanish 


vanish 


vanish 


typical gluon momentum 


asTTlt 







Table 1: Vanishment or non- vanishment of final-state interaction diagrams in the various 
inclusive cross sections. The plus sign "+" shows there is a non-canceled contribution from the 
diagram to the corresponding cross section, while "vanish" shows the cancellation of contribu- 
tions from the diagram to the corresponding inclusive cross section. For details, see subsection 
4.b. 



of the mis-assignment was also studied in Ref. in the open-top region (^/s ^ 2mt) using a 



Monte Carlo generator and with more realistic experimental assumptions. It was shown that 
the effects on the top-quark invariant-mass distribution and on the angular distributions are 
substantial at ^/s = 400 GeV (for = 175 GeV) and increase in magnitude and in complexity 
as the cm. energy is raised. Clearly, in our case, we need more detailed studies to see how 
a similar effect may influence our results. For this purpose, studies based on a Monte Carlo 
generator that produces the fully-differential distribution including the full 0{as) corrections 
near threshold would be necessary. (See also Ref. |12[ for analyses on the radiative color flows 
from t, t, b and b close to the ti threshold.) 

In the course of calculating the final-state interaction effects on the lepton energy- angular 
distribution, we found that the final-state interaction between b and b vanishes (diagram (c) 
and (d)) if we add the corresponding real-emission diagram as well as if we integrate over the 



top-quark energy. This fact was already conjectured in Ref. [|l^ on account of an estimate of the 
Coulomb energy between b and b. Nevertheless the same final-state interaction modifies the top- 
quark energy distribution. As a similar phenomenon one may be reminded of the cancellation 
of final-state interactions (including also those between t and b) in the total tt production cross 
section, despite of the modification of the top-quark three-momentum distribution. In fact the 
cancellation of final-state interactions is a general feature known in a wide class of inclusive 
hard scattering cross sections both in QED and QCD.jJ] 

It may be worth summarizing here at which level of inclusiveness the effects of the final- 
state interactions cancel in the various cross sections in our particular process, top quark pair 
production near threshold and its subsequent decay. This is shown in Table 1. Note that there 
are three typical mass scales involved in this process: the top-quark mass m^, the inverse Bohr 
radius agrrit, and the Coulomb energy between the tt pair a^m^. We show in the table the 



T For example, it is found in Refs. [g6|-|29 that the final-state interaction correction to the invariant-mass 
distribution of W changes sign above and below the distribution peak and that the correction vanishes upon 
integration over the invariant mass. For an enlightening discussion on related problems, see also Ref. js^. 
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typical momentum scale of the gluon in each diagram. 

8 Conclusion 

In this paper we considered the differential distribution of from the semi-leptonic decay of the 
top quark, where the parent top quark is produced in e+e" — > ti near threshold. Particularly, we 
have calculated the final-state interaction corrections (rescattering corrections) to the energy- 
angular distribution of leptons in the top quark decay. Also we have explicitly written down 
the /"•" energy-angular distribution da{e~^e~ ^ tt ^ bl'^iA)W~)/(Pptdxid^i including the full 
0{as) = 0{j3) corrections near threshold. 

We presented numerical studies of the various effects of the final-state interaction correc- 
tions. All numerical results can be understood qualitatively from intuitive pictures. Attractive 
forces between t and b and between t and b modify not only the momentum distribution of 
bW^ or bW~ system but also the top-quark polarization vector and the lepton energy- angular 
distribution. 

• The effect of Coulomb-gluon exchange between b and t, when integrated over the bW~ 
phase-space, can be regarded as a correction to the top-quark production process. The 
effect can be incorporated by modifying the top-quark production cross section and the 
top-quark polarization vector. The top-quark momentum distribution is shifted to take a 
smaller average momentum due to the attraction by b. Also since b is emitted preferably 
in the t spin direction, the attraction generates a cos 9te distribution of the top quark as 
well as modifies the top-quark polarization vector. 

• The Coulomb interaction between b and t causes a non-factorizable correction with respect 
to the production and decay processes of the top quark. It generates an energy-angle- 
correlated correction to the lepton distribution. Namely, the distribution is deformed 
in favour of the kinematical configurations "small Ei and emitted in t direction" or "large 
El and emitted in t direction" , which can be understood as originating from the attraction 
of b in the direction of t. 

• Corrections from the gluon exchange between b and b turn out to vanish when the top- 
quark energy is integrated out. 

Without the non-factorizable effect ^, the angular distribution is dependent only on 
the polar angle from the polarization vector of the parent top quark in its rest frame. The 
final-state interaction brings in another direction into the problem, the direction of t, which is 
a completely new feature in comparison to the decays of free polarized top quarks. 

In order to study the decay properties of top quarks near tt threshold, it is desirable 
to extract the part which is specific to the top-quark decay process alone. In the case of 
semi-leptonic decay, we defined a quantity which depends only on the decay distribution of a 
free polarized top quark. The part which depends on cos 9te and cosOu is dropped using the 
transformation of the l'^ energy and angle which leaves the final-state interaction unchanged. 
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Thus, we recover a differential quantity A{xi,ni-7') dependent only on the lepton energy and 
the lepton angle from the parent top-quark polarization vector. 

This quantity will be useful from the theoretical point of view. It can be calculated from 
the decay distribution of free top quarks without including the bound-state effects or the final- 
state interaction corrections that are typical to the threshold region. Therefore, a variety of 
former studies on free top-quark decays may also be applicable in the tt threshold region, where 
highly polarized top quarks are available with the largest cross sections. 

The authors wish to thank M. Jezabek and J.H. Kiihn for enlightening discussion. 
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